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ABSTRACT 

We report on the Swift-RNI detection of a short burst from the direction of the TeV binary LS I +61°303, 
resembling those generally labelled as magnetar-like. We show that it is likely that the short burst was indeed 
originating from LS I +61° 303 (although we cannot totally exclude the improbable presence of a far-away 
line-of-sight magnetar) and that it is a different phenomena with respect to the previously-observed ks-long 
flares from this system. Accepting as a hypothesis that LS I -|-61°303 is the first magnetar detected in a binary 
system, we study which are the implications. We find that a magnetar-composed LS I -|-61°303-system would 
most likely be (i.e., for usual magnetar parameters and mass-loss rate) subject to a flip-flop behavior, from a 
rotational powered regime (in apastron) to a propeller regime (in periastron) along each of the LS I +61° 303 
eccentric orbital motion. We prove that whereas near apastron an inter-wind shock can lead to the normally ob- 
served LS I +61°303 behavior, with TeV emission, the periastron propeller is expected to efficiently accelerate 
particles only to sub-TeV energies. This flip-flop scenario would explain the system's behavior where a recur- 
rent TeV emission was seen appearing near apastron only, the anti-correlation of GeV and TeV emission, and 
the long-term TeV variability (which seems correlated to LS I +61°303's super-orbital period), including the 
appearance of a low TeV-state. Finally, we qualitatively put the multi-wavelength phenomenology in context 
of our proposed model, and make some predictions for further testing. 
Subject headings: X-rays: binaries, X-rays: individual (LS I +61°303), stars: magnetars 



1. INTRODUCTION 

Besides the highly energetic rotational powered pulsars, 
there is only a handful of other classes of Galactic objects 
known to be emitting until GeV - TeV energies. Some are 
High Mass X-ray Binaries (HMXB). 

The first identified TeV binary system was a 3.4 year pe- 
riod binary hosting a 48 ms radio pulsar in an eccentric or- 
bit around a Be star, namely PSR B 1259-63 (Johnston et al. 
1992, Aharonian et al. 2005). The emission from this ob- 
ject is thought to be associated with the radio-pulsar wind and 
its interaction with the radiation field and/or the material sur- 
rounding the Be-star It shows variable radio to TeV emission 
(Johnston et al. 1999, 2005; Chemyakova et al. 2006, 2009; 
Tam et al. 2011, Abdo et al. 2011b). Oflier two TeV emit- 
ting systems, LS I +61=303 and LS 5039, are botii much 
closer binaries, with orbital periods of 26.5 and 3.9 days, re- 
spectively, hosting a very massive star (Be and O types) and 
a compact object, the nature of which is still unknown. They 
are also both variable, from the radio to the TeV energy range, 
often showing their orbital modulation throughout the multi- 
wavelength energy spectrum (e.g., see, Abdo et al. 2009a,b; 
Aharonian et al. 2006; Albert et al. 2008, 2009; Torres et 
al. 2010; Li et al. 2011). Their X-ray emission showed ks- 
timescale flares (Sidoli et al. 2006; Esposito et al. 2007; Rea 
et al. 2010; Li et al. 201 1), and is, in both objects, character- 
ized by an absorbed power-law spectrum. 

Very recently, other HMXBs emitting at high energies 
have been confirmed. On the one hand, we know of 
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HESS J0632+057 (Aharonian et al. 2007; Hinton et al. 
2009). This system again hosts a Be star in orbit with a com- 
pact object of unknown nature. This TeV binary has an orbital 
period of 320 days (Bongiomo et al. 2011), and again shows 
radio and X-ray variability (Falcone et al. 2010; Skilton et 
al. 2009; Rea & Torres 201 1). No GeV emission has been ob- 
served from HESS 10632+057 yet. On the other hand, Ferm/- 
LAT observations of the gamma-ray source IFGL J1018.6- 
5856 revealed the presence of periodic modulation with a pe- 
riod of 16.5 days (Corbet et al. 2011). Optical observations 
found an 06V(f) star, very similar to that of the gamma-ray 
binary LS 5039, which is coincident with a variable X-ray and 
radio source, leading to the conclusion that IFGL J1018.6- 
5856 is new member of the rare gamma-ray binary class, shar- 
ing in this case several similarities with LS 5039. In the latter 
case, no TeV emission has been reported yet. 

But what is the physical nature of these binaries? Two phys- 
ical scenarios have been put forward (see, e.g., Mirabel 2006). 
On the one hand, a compact object rotating around the mas- 
sive companion star drives relativistic jets as result of accre- 
tion. The gamma-ray binaries would thus be microquasars 
(see, e.g, Bosch-Ramon & Khangulyan 2009 for a review). 
On the other hand, the compact object could be a rotationally 
power pulsar which drives a wind (e.g., Maraschi & Treves 
1981). There are two flavors of models in this case. The inter- 
wind models are those where the acceleration of the electron 
population, which is primary to the gamma-rays observed, is 
originated in the shock region resulting from the interaction of 
the pulsar and stellar winds (see, e.g. Dubus 2006). Instead, 
the physical scenario for high-energy photon production in 
intra-winds models is as follows: pairs are injected by the 
pulsar or inner-wind shocks and travel towards the observer, 
producing Inverse Compton photons via up-scattering thermal 
photons from the massive star. Gamma-ray photons can initi- 
ate an IC cascade due to absorption in the same thermal field 
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Fig. 1. — Sivi/f-BAT light curves of the burst detected in the direction of 
LSI+61°303. 



TABLE 1 

Spectral analysis of the Swift-JiAT burst. 



PL BB Brems 

f 2.0 ± 0.3 

fcT(keV) - 7.5tol iStll 

i?(Km) - 0.27l[!;E;^ 

Flux 5±2 4.5 ±0.7 4.5 ± 0.9 

Fluence 1.4 ±0.6 1.4 ±0.2 1.4 ±0.3 

X^/dof 1.29/14 1.07/14 1.22/14 

Note. — Swift-BfiCT spectroscopy in the 15-50 keV energy range. En'ors 
are at a Icr confidence level for a single parameter of interest. Fluxes and flu- 
ences are in the 15-50 keV energy range and in units of 10~** erg cm"'^ 
and 10~*erg cm~^, respectively. The blackbody radius is calculated at in- 
finity and for a distance of 2 kpc (which is the distance to LS I +61° 303; 
Frail &Hjellming (1991)). 

in both cases (see, e.g., Sierpowska-Bartosik & Torres 2008). 
Models for gamma-rays with pulsars as compact objects have 
been recently discussed by Torres (201 1). 

In the last few decades, two classes of pulsars have attracted 
an increasing interest: the Soft Gamma Repeaters (SGRs) and 
the Anomalous X-ray Pulsars (AXPs). They are two peculiar 
groups of neutron stars which stand apart from other known 
classes of X-ray pulsars. In particular, their X-ray luminosi- 
ties are often larger than the values expected from tapping 
their rotational power reservoir (L ^ lO'^^erg s^^), and they 
show no evidence for a companion star which could power 
this strong emission via accretion. Their rotational periods 
(2-12 s) and period derivatives (^ 10^^"^ — lO^^^B), for most 
of the ^20 sources known to date, point to a magnetic field 
of ~ 10^'' — 10^^ G, which is currently believed to be re- 
sponsible for their peculiar emission properties (see Rea & 
Esposito 2011 and Mereghetti 2008 for recent reviews). At 
present, in fact, the model which is most successful in ex- 
plaining SGRs and AXPs emission is the "magnetar" model: 
these objects are thought to be strongly-magnetized (isolated) 
neutron stars emitting across all wavelengths via the decay 
and the instabilities of their high B-fields (Duncan & Thomp- 
son 1992; Thompson & Duncan 1993). 



The most peculiar and intriguing property of SGRs and 
AXPs are their outbursts and flares, which are at variance 
with any other bursting event observed so far in other com- 
pact objects. In particular, the unpredictable flaring activity 
of magnetars can be phenomenologically classified in a few 
types: 

1. X/'y-ray short bursts: these are the most common and 
less energetic flaring events, they have short duration 
(^0.01-0.2 s), thermal spectra, and peak luminosities 
of ^ 10"^^ — lO"*^ erg s^^, and they can occur randomly 
as single events or in a bunch (see recent examples in 
Rea et al. 2009; Kumar et al. 2010). 

2. Intermediate flares: they are intermediate both in dura- 
tion and luminosity between short bursts (1) and the Gi- 
ant Flares we discuss below. They have durations rang- 
ing between ^2 - 60 s and luminosity of ~ Iff*^ — 10 
erg s^^. Sometimes intermediate flares last longer than 
the pulsars spin periods and shows a clear modulation at 
the star rotational period (see recent examples in Israel 
et al. 2008; Mereghetti et al. 2009). 

3. Giant flares: they are by far the most energetic (^ 
10*^ — 10^^ erg s~^) Galactic events after supernova 
explosions. We have detected so far only three of these 
events, which are characterised by a very luminous hard 
peak lasting less than a second, and by decaying tails 
lasting 100-500 s where the spin period of the neutron 
star is clearly visible. 

4. Outbursts: they are enhancements of the multi-band 
emission of SGRs and AXPs by a factor of 5 - 1000, 
with a typical total energy release of ^ 10""' — 10''^ erg. 
The increase flux level of the source may last from a 
few months up to several years (see Israel & Dall'Osso 
201 1 and Rea & Esposito 201 1 for a recent review on 
the flaring activity of magnetars). 

All of the above flaring events are peculiar to, and defining 
the magnetar class. 

In this paper, we first provide a complete analysis of a 
magnetar-like short burst detected by the Burst Alert Tele- 
scope (BAT; Barthelmy et al. 2005) onboard Swift from 
the direction of the TeV binai-y LS I +61° 303. The event 
prompted several instant notices (Astronomer Telegrams and 
Gamma-ray bursts Coordinate Network circulars): the Swift- 
BAT trigger was reported by De Pasquale et al. (2008), the 
5w//f-BAT localisation by Barthelmy et al. (2008), and Evans 
et al. (2008) noticing that a quasi-simultaneous follow up 
with the Swift-XRT unveiled only one source within the burst 
error circle: LS I +61°303. Furthermore, brief interpreta- 
tional comments and reports of multi-wavelength observa- 
tions already discussing the magnetar-like features of the burst 
were posted (Dubus & Giebels 2008, Rea & Torres 2008, and 
Munoz-Arjonifla et al. 2008). 

To try understanding the nature of this short burst, the 
Swift-BALT data were re-analysed, as well as those of 
various RXTE and Chandra observations of the field of 
LS I +61° 303. In particular, the data coflected by RXTE si- 
multaneously with its Proportional Counter Array (PCA) and 
High-Energy X-ray Timing Experiment (HEXTE; the analy- 
sis of these data is presented here for the first time) allows us 
to explore a possible alternative origin for the burst, namely, 
the possibiHty of it being the result of a spectral evolution 
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TABLE 2 

RXTE-UEXTE COUNTERPARTS TO THE SOFT X-RAY FLARES. 



Fig. 2. — S»'iyf-BAT positional error circles over-imposed to the Chandra 
image of the field (smoothed with a Gaussian function with a FWHM of 
3"). Cyan, green and yellow circles report on the best position derived by 
Barthelmy et al. (2008), and our analysis using the T90 and the T(tot) of the 
burst (see text for details), respe ctive ly. The source within the red circle is 
the Chandra source discussed in 92.31 



in a longer (normal) flare. All of our observational analysis 
are presented in S|2] In f|3l we compare the 5w//f-BAT burst 
with those observed fromknown magnetars, showing in de- 
tail the striking similarity with the detection of interest. In f|4] 
we confront the hypothesis that LS I +61° 303 might indeed 
host a magnetar, and start exploring its implications. In f|6] we 
discuss our results and the proposed model resulting from the 
previous sections on the light of the multi-band variable emis- 
sion of LS I +61° 303 reinterpreting all observations from ra- 
dio to TeV gamma-rays when we place ourselves under the 
possibility of having at hand evidence for the first magnetar 
in a binary system. 

2. OBSERVATIONS AND ANALYSIS 
2.1. Swift- BAT: a magnetar-like burst 

The Burst Alert Telescope (BAT; Barthelmy et al. 2005), 
the hard X-ray detector on board Swift (Gehrels et al. 2004), is 
a highly sensitive coded mask instrument optimized in the 15- 
150 keV energy range. It was specifically designed to catch 
and study prompt emission from gamma-ray bursts and other 
interesting high-energy transients. 

On 2008 September 10 at 12:52:21 UT, the Swift-BAT 
triggered on a burst in the direction of the gamma-ray bi- 
nai-y LS I +61°303 (trigger #324362; De Pasquale et al. 
2008; Bai-thelmy et al. 2008), see Fig. [T] The data cal- 
ibration and reduction was performed using the standard 
BAT analysis software distributed within FTOOLS under the 
HEASOFT package (version 6.9), and the latest CALDB re- 
lease (2011-03-03). The mask-tagged (i.e. background- 
subtracted) counts of the source were extracted from the de- 
tector pixels illuminated by the source by using the mask- 
weighting technique. The mask-weighting factors were cal- 
culated by BATMASKWTEVT using the gi-ound-calculated po- 
sition of Barthelmy et al. (2008, GCN #8215; RA = 40.101°, 
Decl. = 61.210°, about 2.9 arcmin form the BAT on-board 
position reported by De Pasquale et al. 2008, GCN #8209; 
see also Fig.|2]l. Fig. |3]present the 5w//f-BAT image of the 
burst, together with images of the same region obtained with 
different instruments, as discussed below. 

Mask-tagged light curves were created in the standard 4 en- 
ergy bands, 15-25, 25-50, 50-100, 100-350 keV (Fig.[T]) at 
64 ms time resolution. The burst is visible in the first two 
bands (15-25 and 25-50 keV), while no significant excess 
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Note. — The average count rate, significance and reduced in the 15- 
250 keV energy band of the HEXTE data during the ks-timescale flares. 

is observed above 50 keV. The total duration of the event in 
the 15-50 keV band is Ttot ^ 0.31 s, while the Tgo dura- 
tion is 0.24 ± 0.05 s. These values were computed by the 
BATTBLOCKS task (based on the Bayesian Block algorithm; 
Scargle 1998) from a light-curve with 1 ms bin size. 

We extracted a 15-50 keV sky image and performed over 
the Tgo duration of the burst a blind source detection with the 
tool BATCELLDETECT. This script performs a least-square fit 
of peaks in the map to the BAT point-spread function (a two- 
dimensional Gaussian) using the local rms noise to weight 
the pixels in the input map. A single, highly-significant 
(11. Oct) point-like source was detected at the best-fit coor- 
dinates RA = 40.1119°, Decl. = +61.2322° (essentially 
identical results were obtained using the total duration of the 
event: RA = 40.0962°, Decl. = +61.2362°; signal-to-noise 
ratio 11. Oct; see Fig.l2]i. 

The Point SpreaciFunction fit using batcelldetect 
also yields a formal uncertainty based on the least-square co- 
variance matrix. At Icr this was 1.06 arcmin for the run using 
Tgo (and 1.04 arcmin with Ttot)- However, because neigh- 
boring pixels in the coded mask images are inherently corre- 
lated, this error is known to be a poor estimator of the true 
uncertainty for high signal-to-noise detections. Therefore we 
adopt a more conservative figure of 1.4 arcmin at la (90%: 
2.1 arcmin) following the prescriptions of the BAT calibra- 
tion reportsjand including also a 0.25-arcmin systematic er- 
ror (see Tueller et al. 2010). This position is consistent with 
that reported by Barthelmy et al. (2008) (albeit the center of 
the uncertainty visibly moves) and with that of LS I +61° 303 
(the angular separation being ^0.6 arcmin); see Fig. |2] No 
source is detected in the BAT image excluding the burst mter- 
val. The burst spectrum was extracted over the Ttot interval 
(we extracted it also from Tgo finding consistent results). The 
results of the spectral analysis are summarized in Table [T] To 
investigate the spectral evolution as a function of time, we 
computed a hardness ratio, but no spectral variations in the 
hardness ratio were visible. 

2.2. Analysis of RXTE-HEXTE observations: can the 
Sv/ift-BAT short burst be the result of spectral 
evolution ? 

The RXTE-PCA observation nearest to the ^wZ/f-BAT burst 
is located about 6 hr after the burst, with a total exposure 
time of 1492 s. We have analyzed this observation before, 
and no unusual behavior (and no flare) was found (see Ray et 
al. 2008, Torres et al. 2010). 

In order to further investigate a possible connection be- 
tween the short magnetar-like burst seen by Swift-BAJi (see 
j j2.1[ l, and the typical ks-timescale bursts often observed in 
LST^ +61° 303 (Sidoli et al. 2006; Rea et al. 2010; Smith 



^ See http://swift.gsfc.nasa.gov/docs/swift/analysis/bat_digest.htmI 
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Fig. 3. — From left to right: Swift-BKT image of the short burst (see j ]2.1| , superimposing the spatial accuracy of our position determination (1.4'); Swift-XKT 
image of all data collected so far on LS I +61° 303 (165 ks), with the 1 A Swift-BKT error circle superimposed); Chandra ACIS-I 50 ks image of the field of 
LS I +61° 303 with the short burst positional accuracy over-imposed. 

numbering scheme) in the second cluster (B) started to func- 
tion abnormally since 1996 March 6. Because of this, only 
Detector 0, 1, 3 in cluster B are used in our analysis till (in- 
cluding) 2009 December 13. After that date, and as suggested 
by the RXTE team, we take data from cluster A as the source 
and data from cluster B as the background measurements. We 
select time intervals where the source elevation is > 10° and 
the pointing offset is < 0.02°. HEXTE lightcurves were gen- 
erated, deadtime corrected, and background subtracted using 
REX script. 

In the bottom panel of Fig. |4] we show the 15-250 keV 
RXTE-HEXTE light-curve of LS I +61° 303 in 16 seconds 
time bins. For comparison we show in the top panel of 
the same figure also the 3-lOkeV RXTE-PCA light-curve of 
LS I +61° 303 in 16 seconds time bins. The latter enhances, 
because of the addition of newer data, the results presented 
by Li et al. (201 1). In particular, it can be seen there that we 
discovered a sixth flare out of the PCA monitoringr] 

Since, as we explained above, the overall HEXTE 
lightcurve is produced with different instruments at differ- 
ent times -cluster B as source and background before (and 
including) 2009 December 13 (MJD 55178), and cluster A 
as source with cluster B as background after that date- the 
average of the lightcurve changes, which is obvious to the 
eye in the bottom panel of Fig. [4] Fitting a constant to the 
HEXTE lightcurve before MJD 55178 yields an average count 
rate of 0.11 ± 0.03, and a reduced = 1-29 (12915 d.o.f). 
Data accumulated for more than two years (2007 Septem- 
ber - 2009 December) generates a detection significance for 
LS I +61°303 of 4.4cr in the 15-250 keV band. The aver- 
age count rate after MJD 55178 is 4.99 ± 0.03 and a fit to a 
constant yields a reduced = 1-76 (10745 d.o.f). 

To investigate whether there are possible flares in RXTE- 
HEXTE, we looked at data points with significance above 4(t 
over the average count-rate. In the 15-250 keV band, 18 out 
of 23662 bins points have a significance large than Aa. None 
of these corresponds to any of the flares detected by the PCA 
in the soft X-ray band, and we associate them to statistical 
fluctuations. A plot of the number of data points with a given 

* The data set additionally analyzed in this report covers the time span 
from 2010 September 5 to 2011 February 27 (MJD 55444^55619), and covers 
nearly another half year of data recently released on the HEASARC website. 
It includes 50 RXTE pointed observations identified by proposal numbers 
95102 and 96102, providing a total exposure time of 70 ks on the source 
beyond that reported by Li et al. (201 1). We follow exactly the same analysis 
chain as reported by Li et al. (201 1) in analyzing the newest data. The general 
properties of the new, sixth flare we discover here are very similar to those 
of the fourth flare reported by Li et al. (201 1). Its power spectrum shows no 
evidence for the existence of any structure. 
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Fig. 4.— Top panel: RXTE-PCA long-term light-curve of LS I +61°303 
in 3-10 keV. Bottom panel: RXTE-UEXTE light-curve in 15-250 keV. The 
highlighted regions are the 5 flares found in PCA data (Li et al. 2011) plus 
a sixth additional one we uncovered in this paper as a result of continuous 
analysis of the new data. 

et al. 2009; Li et al. 2011) we searched for a possible hard 
X-ray counterpart to the ks-timescale flares observed by the 
RXTE-PCA, taking advantage of the simultaneous observa- 
tions performed by the RXTE-HEXTE. The main idea was to 
search for any short/hard X-ray burst in coincidence with the 
longer flares observed in the 3-10 keV energy range (Li et al. 
2010), assuming that a strong spectral softening might occur 
the first seconds after the flare emission. 

RXTE-HEXTE consists of 8 detectors with a total area of 
about 1600 cm^. The 8 detectors are split up in two clusters 
of 4 detectors each. During normal operations, each cluster 
is alternately pointed on and off source, generally every 16 
or 32 s, to provide a nearly real-time background measure- 
ment. We use here HEXTE data covering the period between 
2007 September and 201 1 February. The data encompass 418 
HEXTE pointed observations, providing a total exposure time 
of 621 kson LSI+61°303. 

In the lightcurve analysis of HEXTE we used "Standard 
Modes" data with 16 seconds accumulation time, and built 
the lightcurve in the 15-250 keV range. Data reduction was 
performed using the HEASoft tools, and the data were fil- 
tered using standard HEXTE criteria. At the end of 2006, the 
first cluster (A) of detectors was fixed to always measure the 
source of interest and no background measurements were car- 
ried out since then. From 2009 December 13, cluster B was 
fixed to always measure the background and no source detec- 
tions are available with it. Additionally, Detector-2 (in the 0-3 
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Fig. 5. — Lightcurve of the 6 flares in 3-10 keV (PCA) as observed at 15-250 keV by HEXTE. All lightcurves are binned at 16 s. 



count rate would show no deviation from a Gaussian fitting. 
Moreover, we divided the 15-250 keV energy band into 15- 
60 keV and 60-250 keV and found that all points above 4(t 
significance in the larger energy range, 15-250 keV, are no 
longer that significant in 15-60 keV and/or 60-250 keV. FigJS] 
shows, corresponding to the 6 PCA flares marked in Fig. BT 
a zoom of the HEXTE lightcurve. The average count rate, 
significance, and reduced of a constant fit to each flare 
are listed in Table |2] and show that all 6 PCA flares are not 
significantly detected in the 15-250 keV band. 

If any of the PCA flares would have been preceded by a sim- 
ilar burst to the one detected by 5w//f-BAT, HEXTE would 
have seen them. To show this, we have simulated a burst 
observed by HEXTE with ^wZ/f-BAT parameters. Using the 
HEXTE response, we found a significance of 4.64a for it. If 
the Swift-BAI burst would have been the result of spectral 
evolution, and related to the longer timescales flares usually 
detected from LS I +61° 303, then HEXTE would have spot- 
ted flares at least in any of the six instances in which we had 
PCA and HEXTE simultaneous coverage. This has not hap- 
pened, HEXTE have not seen any of the PCA ks-timescale 
flares, and we conclude that their nature is different from the 
^wi/f-BAT detection. 

2.3. Chandra ACIS-I: search for a serendipitous magnetar in 
the field of LS I +61° 303 

We analysed a ^50 ks observation of LS I +61° 303 
performed with the Advanced CCD Imaging Spectrometer 
(ACIS) instrument (ObsID 10042) starting on 2006 April 07 
22:08:59 (UT), in VERY FAINT (VF) timed exposure imag- 
ing mode (see also Paredes et al. 2007). The source was 
positioned on the back-illuminated ACIS-I3 CCD. Standard 
processing of the data was performed by the Chandra X-ray 
Center to Level 1 and Level 2 (processing software DS 8.0). 
The data were reprocessed using the CIAO software (ver- 
sion 4. 1 .2). We used the latest ACIS gain map, and applied 
the time-dependent gain and charge transfer inefficiency cor- 
rections. The data were then filtered for bad event grades 
and only good time intervals were used. No high back- 
ground events were detected. The final net exposure time was 
49.105ks. 

We used the CIAO celldetect tool to search for sources 
by summing counts in square cells in the dataset, and compar- 



ing the counts to those of "background" cells. At each point 
where a cell is placed, a signal-to-noise ratio of source counts 
to background counts is computed. We placed a detection 
limit of S/N=4, and found only one source marginally com- 
patible with the refined position of the 5w//f-BAT burst (see 
j pl) at RA = 40.095081°, Decl. = +61.258468° (with a 
Icr error of 1.5"on the position), with a total of 0.3-lOkev 
counts of 50 ± 9 counts (see also Fig. [2]), translating in a 
count-rate of 0.0010 ± 0.0002 counts s^ VNo radio, infrared 
or optical counterparts have been detected for this source de- 
spite the deep archival observations covering this field of view 
(see Munoz-Arjonilla et al. 2009; this source corresponds to 
their #12 of Table 3). Thus, formally we can not exclude that 
the faint X-ray source detected by Chandra at the limit of 
the 1(7 positional uncertainty of the burst (see Fig. |2]l might 
be the magnetar responsible for the short burst observed by 
^wZ/f-BAT, hence independent from LS I +61° 303: Assum- 
ing a thermal spectrum of ^ 0.3 keV, typical of a magnetar 
in quiescence (see Rea & Esposito 2011 for a review), and 
an absorption column density of 9 x lO^^cm^-^ (relative to 
the whole Galactic value in the direction of the source, from 
the HI maps from Dickey & Lockman (1990)) we derived a 
0.3-lOkeV observed flux of ^ 6.1 x lO^^^erg s^^cm^^. As- 
suming the source is located at the end of the Milky Way, 
at lOkpc distance, the corresponding luminosity would be 
^ 2.7 X lO'^^ergs^^, consistent with it being a magnetar in 
quiescence. Given that the number of TeV binaries is a hand- 
ful, and the number of magnetars in the Galaxy that have been 
detected by our experiments is also a handful, the probability 
that both are seen at ^ 1.4' from each other seems a priori 
low. A precise number can not be computed without fur- 
ther assumptions at many levels (population distribution, total 
number of sources, etc) which would probably make mean- 
ingless the result since we can never rule out a single random 
coincidence. The possibility is thus left open, although it does 
not seem to be one that would be reasonably preferred. 

3. THE SWIFT-BAT BURST IN THE CONTEXT OF MAGNETAR 
EMISSION 

The properties of the burst observed by ^wZ/f-BAT (a very 
short duration an d a thermal spectrum with a temperature of 
~ 7.5 keV; see S2.1i are typical of magnetars (see Aptekar 
et al. 2001, WooHs & Thompson 2006; Mereghetti 2008), 
and at variance with other kinds of frequently observed flares 
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as Type I bursts from X-ray binaries (which last ^100 times 
longer) or Gamma-ray bursts (which have harder spectra and 
are much more energetic). In particular, the burst flux (see 
TableriJ at a distance of 2.3 kpc (as for LS I +61°303) implies 
al5-50keV luminosity of^ 2 x lO'^^erg s^^. The luminosity 
of this burst is in the lower end of the distribution of short 
bursts from magnetars, in line with the bursts observed in the 
AXPs (see Gavriil & Kaspi 2002; Woods et al. 2004) which 
are usually slightly less powerful (and lasts longer) than the 
bursts observed from the canonical SGRs (Gogus et al. 1999, 
2001; Israel etal. 2008). 

The 'relatively' low intensity of the burst, and the single 
burst that has been found in decades of observations towards 
LS I +61° 303 and in X-ray all sky surveys, perhaps suggest 
that the magnetic neutron star producing it has a magnetic 
field at the lower end of the typical magnetar regimes (hence 
around - 5 x 10" G as the case of PSRB1846-0258; Gavriil 
et al. 2008; Kumar & Safi-Harb 2008), or it is a rather old 
magnetar (Myrs; see e.g. Perna & Pons (2011)). Note that in 
the first scenario a high rotational power can be present (E ^ 
Iq35 _ lo^^ergs^^), while in the latter scenario the object 
would be a rather slow pulsar with little energy stored in its 
rotation (see e.g. Rea et al. 2010 for one of such examples). 

However, in both cases the X-ray emission of LS I +61° 303 
can easily be not magnetar-like (e.g. due to resonant cy- 
clotron scattering of a hot surface through a very dense mag- 
netosphere; Thompson, Lyutikov & KuUcarni 2002). In fact, 
young magnetars with a relatively low B field and high ro- 
tational power have X-ray spectra dominated by non-thermal 
processes due to particle acceleration and shocks from their 
strong winds, while old magnetars in quiescence are rather 
faint (10'"^ — lO'^^erg s^^), having dissipated most of their 
magnetic energy, and emit mostly thermally from their sur- 
face. Thus, from the above considerations and what we will 
explain in detail in the following sections, we believe that in 
LS I +61° 303 we might be witnessing a high-ii^ magnetar 
with a magnetic field of the order of 10^'' — 10^* G. 

4. MAGNETAR IN A BINARY SYSTEM 

If we entertain the hypothesis that the origin of the 
magnetar-like event reported here could be LS I +61°303, 
an effort to understand what are the consequences of the exis- 
tence of a magnetar in this eccentric binary system is in order. 
How will the observed multi-wavelength phenomenology be 
generated in such a case? The rest of this paper is devoted to 
analyze these questions. 

4.1. Physical radii 

To start considering these issues, we shall introduce several 
radii (measured from the neutron star), which will represent 
the relative strength of the system's components and compare 
those with the position of the light cylinder, Ric (see, e.g., 
lUarionov & Sunyaev 1975; Davies & Pringle 1981, Lipunov 
etal. 1994). The latter is 

Ri,= ^o^ 4.77 X 10^ (P\ cm, (1) 
2n V Is/ 

where P is the period of the pulsar, and we have adopted for it 
a scale of 1 s, typical of the observed magnetars. The Alfven 
or magnetic radius, Rm, will be defined as the distance in 
which the magnetic field starts to dominate the dynamics of 
the in-falling matter. Thus it can be implicitly defined by 



the equality between the energy densities, which is attained 

at R,n, 

1 , 

Here, Vf is the free fall velocity of the accreting matter onto 
the neutron star of mass Mns, 



2GMns 



R 



(3) 



B is the magnetic field in the inner magnetosphere, assumed 
to be of dipole type with neutron star radius Rns 



B{R) - B„ 



1i 



(4) 



where R is measured from the neutron star center, and p is 
the density of the accreting matter The value of the latter is 
defined as 

where Mace denotes the rate at which matter is accreted, say at 
a distance r from the companion star, and can in first approxi- 
mation be obtained using the Bondi-Hoyle-Littleton approach 
as (see, e.g., Bondi & Hoyle 1944) 



Maec(r) = 

4 \ r 



(6) 



In the latter formula, Af* denotes the stellar mass-loss rate 
(and again r is measured from the companion, not the neu- 
tron star). Finally, Reap is the gravitational capture radius for 
a neutron star of mass Mns, as determined by the relative ve- 
locity of the neutron star with respect to the matter, Vrei, using 
the orbital and the stellar- wind velocity 



cap 



•2GMns 

rel 



(7) 



4.2. Regimes under the influence of a polar wind 

The stellar wind velocity description will first be assumed 
having the typical form of a radiatively-driven outflow from a 
high-mass star (e.g.. Castor & Lamers 1979), 



K„ - + {Voo - K)) 1 



R, 



1 



R, 



r J \ r 

(8) 

where is the stellar radius, Vq ^ 0.01^00, /3 ~ 1, and 
Foo ^ 1000 km s-^ (Lamers & Casinefli 1999). Note that in 
Eq. (ISll, r is measured from the companion, thus, at R,n, we 
will have r{Rm) = d — with d being the system separa- 
tion. In an elliptic orbit of eccentricity e, the system separa- 
tion is given by 



a(l — ecos(e)), 



(9) 



where a is the semi-major axis and e is the eccentric anomaly 
(see, e.g., Hilditch 2001). Considering that d ^ Rm, we 
neglect R„i in favor of d in the definition of r{R„i) for the 
subsequent calculations. We have checked, however, that 
not assuming this approximation would correct the results we 
achieve in a few percent while significantly complicating the 
algebra. We also consider that the polar wind terminal ve- 
locity, Voo, dominates over the orbital speed (this is explicitly 
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shown below). Taking this into account, the capture radius is 
given by 



Rrn.T) 



2GM„ 

^2 



= 3.73 X 10 
1 - 0.69 



10 / ^-^ns \ I ^nc 



IAMqJ \10^ cm s-i 



WRqJ (i012 cm) 
(1 — ecos(e)) 1 cm. 



(10) 



We shall consider that a pulsar acts normally, i.e., being ro- 
tational powered and driving out a relativistic wind, if at the 
capture radius the pressure of the pulsar wind (cosmic rays 
and magnetic field), given in terms of its spin-down luminos- 
ity Lsd as 

exceeds the pressure of the stellar wind gas behind the shock 
(Pw — Pwyw)^ s^^' ^-S-' Illarionov & Sunyaev (1975). In this 
case, the matter is swept away beyond Reap and the system is 
so-called ejector (Lipunov et al. 1994). TTie pressure condi- 
tion can be written as 



Lsd> 'iTT RlapCPwV^. 



(12) 



This entails a condition onto the pulsar period (for other mag- 
nitudes fixed): in order for the pulsar to stop acting as an ejec- 
tor, the period should be larger than 



Rout- 



ejector 



1 S 



>3.09 



B ( M. 



1014 G 



g 



10^ cm 



3/4 



10^^ cm 



-1/4 



1/2 



lO" cm 



3/2 



I AM, 







-1/2 



X (l-ecos(e))i/2 
X ( 1 - 0.69 ^* 



(1 — ecos(e)) 



l\3/J/4 



(13) 



In these cases, with decreasing radius (i? < Reap), the pres- 
sure of the accreting matter grows as i?-^/^ whereas the Ppsr 
does it as R^^. Because of this, matter fall from Reap be- 
gins, penetrating within the light cylinder up to the point (if 
any) in which the quick rise of the dipolar magnetic field in 
the magnetosphere produces a pressure (oc R^^) able to stop 
it (i.e., at the Alfven radius R„i). The pulsar no longer has 
a magnetosphere and cannot generate a relativistic wind any 
longer, leaving the ejector phase. In order to move out of this 
stage, and to ignite the normal (rotationally-powered) pulsar 
again, an unscathed magnetosphere should be recovered, and 
thus the condition for it to happen is that the period be smaller 
than that needed to have Rm ~ Rie- 
The value of Rm can be obtained using previous expres- 



Eq. ; 



i?„ =2.1 X 10' 



B 



4/7 



AL 



-2/7 



10" G/ \10i8gs- 



8/7 



4/7 



lO^cms-1/ V 10^2 cm/ 

12/7 / ,,f X -5/7 



( -^^) ( (1 - e cos(6))4A 
Vl06 cm J \1AMqJ ^ ^ " 

f R* 



X 1 - 0.69 



V lOi?0 J V 1012 cm 

_l\8/3/7 

(1 — ecos(e)) 1 



cm. (14) 



Note that whereas Rie depends (linearly) on the neutron 
star spin period, does not depend on P at all. The line 
Rm = Ric thus entails a condition onto P. Neutron stars of 
sufficiently small periods have an unscathed magnetosphere 
even when in eccentric orbits with high mass-loss rate stellar 
companions, and thus behave as a normal pulsar. In general, 
in order to re-ignite the pulsar the condition over the period 
(considering all other magnitudes fixed) reads 



Pinto— 



ejeetor 



1 S 



<4.5 



B 



10" G 



4/7 



-2/7 



IQis gs-1 



( )"'( 

VlOScms-iy V 



4/7 



1012 cm/ 



/ Rn 



12/7 



-5/7 



VlO'^cmy V 1-4*^0 

X (l-ecos(e))4/7 



X 1 - 0.69 



VlQis 



VlOi?0 J ViU^^ cm 

l\ 8/3/7 



-l\ 

(1 — ecos(e)) 1 



(15) 



Using the fiducial values in Eqs. ( [T3] l and ( [T5] l we see that 
a neutron star with an spin period usually measured for mag- 
netars would be right at the transition range. Small changes 
in the mass-loss rate, for instance, can make the neutron star 
flip-flop from accreting within the magnetosphere to behaving 
as a rotational powered pulsar. 

To simplify the reasoning that follows, we will consider 
that the condition Rm = Rie not only establishes the into- 
ejector condition over the period P, but also the out-of-ejector 
regime. Indeed, the out-of-ejector condition obtained with Eq. 
( 13 I is easier to fulfill (i.e., slightly smaller periods can fulfill 
itTTor equal values of the other magnitudes involved) than that 



obtained with Eq. ( 15 1, whereas the shape of the constraint is 
similar. Fig. [6]compares these two constraints over P, differ- 
ing by a factor of ^ 1.57 one from another 

4.2.1. Orbital eccentricity ejfects 
under the influence of a polar wind 

Fig. It] shows the condition Rm — Ric as given by Eq. 
15 I for oifferent values of the neutron star (surface) magnetic 



field. We adopt one order of magnitude up and one down for 
the possible variation in these parameters out of their fiducial 
values noted in Eq. ( [T5] l. We assume that the radius of the 
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M (1 0" g s ') Eccentric anomaly (c / ° ) 



Fig. 6. — Comparison of the out of ejector conditions written in E g. )13) - 
in violet- and the constraint of the equality Rm = Ric, given by Eq. \\~t\ -in 
red. The latter is always a more restrictive condition onto P, for equal values 
of the other magnitudes involved. Two cases are shown for surface magnetic 
fields of 5 xlO^'* G and 5 xlO^^ G. The eccentricity is assumed equal to 
zero and the LS I +61° 303 semi-major axis is adopted in this example. 



100 




1.0 

{M./10"gs') 



10.0 



Fig. 7. — Phase space conditions for accretion regimes around a highly 
magnetic neutron star under the influence of a polar wind of terminal velocity 
Voo and mass-loss rate A/» . It is assumed that the radius of the neutron star is 
1 0^ cm and its mass is 1 .4 M0 . An example of the phase space separation is 
marked for the case of surface magnetic field of 5 X 10^'' G. The eccentricity 
is assumed equal to zero and the LS I +61°303 semi-major axis is adopted 
in this example. 

neutron star is 10^ cm and its mass is 1.4 M©. We adopt 
a semi-major axis value of 6 x 10^^ cm and a radius of the 
massive star equal to 10 R©, consistent with measurements 
of the LS 1 +61°303 system (see e.g., Casares et al. 2005, 
Grundstrom et al. 2007). 

An example of the different physical behaviors of the sys- 
tem is noted in Fig. [vjby the Hne with _B = 5 x lO^'^ G, 
separating the region of the plane, for that specific value of 
magnetic field, where the accretion is halted within the mag- 
netosphere (in light green, where Rm < Ric), from that in 



Fig. 8. — Ratio of the _Rm-values along an eccentric orbit (under the influ- 
ence of a polai' wind of terminal velocity Voo ) with respect to that attained at 
periastron, for different eccentricity, as a function of the eccentric anomaly. 

which the systems acts as a rotational powered system (in 
light yellow, where > Ric). This separation is valid at 
each of the possible (neutron star surface, dipolar) magnetic 
field strengths, which form a continuum throughout the plot. 

Fig. |7] and Eq. ( 15 1 show, as mentioned in the last Sec- 
tion, that the observea values for periods and magnetic fields 
of known magnetars, together with fiducial (and commonly 
adopted) values of the stellar (polar) wind velocity and mass- 
loss rate of LS 1 +61° 303 would make, assuming a circular 
orbit, Rm ^ Ric- For example, a system hosting a pulsar with 
P ^ 7 s, where measured periods of magnetars cluster, and 
the fiducial values for the properties of the wind of the mas- 
sive star, Vw — 10* cm s^^ and A/* ^ 10^* g s^^, would be 
right at the line representing _R,„ = Ric for a surface magnetic 
field B = 5 X 1013 G. 

However, one has to take into account that the orbit of 
LS 1 +61° 303 is not circular, and its eccentricity have been 
quoted in the range of 0.55 to 0.72 (again see, e.g., Casares 
et al. 2005, Grundstrom et al. 2007, Aragona et al. 2009). 
Rm then becomes a function of the orbital position, and can 
be represented, for a given value of eccentricity, as a func- 
tion of the system's eccentric anomaly. Fig. |8]shows that the 
ratio of the i?.m-values along an eccentric orbit with that at- 
tained by Rm at periastron. For the quoted eccentricity, Rm 
is a factor between 2 and 4 smaller at periastron than what it 
is at apastron. Then, if the magnetar-composed system fulfills 
the condition Rm > Ric at apastron and the inner magneto- 
sphere thus behaves as that of a rotational powered pulsar, it 
will likely fulfill Rm < Ric at periastron. This is graphically 
shown in Fig. |9] which shows, for e ~ 0.6, the region of 
the phase space in which a magnetar hosted in LS 1 +61° 303 
would change behavior. Fiducial values of all parameters, un- 
der the assumption of a dominant polar wind, would position 
the system right in the middle of the flip-flop regime, with the 
putative magnetar in LS 1 +61° 303 acting as a rotational pow- 
ered system in apastron {Rm is larger than the light cylinder 
in apastron), and accreting within the magnetosphere in pe- 
riastron {Rm is smaller than the light cylinder in periastron). 
For usual magnetar parameters, this is the natural solution for 
the LS 1 +61° 303 behavior along its orbit. 
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4.2.2. Kind of accretion under the influence of a polar wind 

In order to consider further the possible accretion scenario, 
we will define two additional radii. Let us consider first the 
rotational velocity of the magnetosphere at the position 
which is given by 

Vrot = p , (16) 

and the Keplerian velocity of the infalling matter, which is 
instead given by 



100 



Vi 



Kep 



Rr, 



(17) 



These two velocities allow for the definition of the co-rotation 
radius, Rcor, when Vxep = Vrot, which can be written as 



Rr.nr 



fGMnsP^ 

\ 47r2 
1.67 X 10^ 



1/3 



system always accreting 




system always rotationally powered 



0.1 



1.0 

(M./10'"gs 



10.0 



1.4M 



o 



1/3 



P 

Is 



2/3 



cm. (18) 



The co-rotation radius represent the position of the centrifugal 
barrier for the in-falling material created by the neutron star 
rotation. i?,„ > Rcor when 



P 

Y7s 



< 1456.3 



B 
10" G 




-3/7 



10^ cm 



V V 1012 cm/ V 



R„ 



10^ cm 



18/7 



Fig. 9. — Examples of the region of the phase space for a system similar 
to LS I +61°303 where the system would change behavior from accreting 
within the magnetosphere (in periastron) to be a rotationally powered pulsar 
(in apastron). Each of these phase space regions producing this flip-flopping 
behavior are located within the two lines corresponding to equal magnetic 
field label; for which several pairs are shown. An example with B = 5 X 
10^^ G is explicitly depicted (the central white region of the plot is the flip- 
flopping area; the top green part corresponds to always-accreting systems, 
and the bottom light-yellow one to neutron stars always acting as a rotational 
power pulsar along the orbit). The flip-flopping region moves up in the plot 
for higher magnetic fields (an example with B = 10^* G is shown with red 
lines). The neutron star is assumed to be subject to the influence of a polar 
wind of terminal velocity Voo ■ 



1.4M0 
1 - 0.69 



-37/14 



(l-ecos(e))^/^ 



(1 — ecos(e)) 



-2/3 



(20) 



R* 



WRq 

(1 — ecos(e))^"^^ 



1012 cm^ 

12,3/7 



In general though. Reap > Rm implies, e.g., a condition upon 
the stellar wind velocity as a function of the other parameters 
in phase space. 



(19) 



That is, unless values far-off from the fiducial ones are in- 



10* cm s^i 



< 1.2 



B 



-2/11 



M, 



1/11 



1014 cy \ioi8gs-i 

2/11 / \ -6/11 



1012 cm/ V 106 cm 

( ^^V^"(l-ecos(e))-2/ii 



when Reap < 
when Reap > 



X 1-0.69 



lORi, 



1012 cm 



voked, the inequality in Eq. ( 19 1 is always respected. Particu- 
larly, when the period is suchthat Rm is less, but of the same 
order than Ric, the system is halting accretion at distances 
from the neutron star which always results in super-Keplerian 
velocities. 

When the accretion is halted in the magnetosphere, i.e., 
Rm < Ric, at a position which is super-Keplerian, Rm > 
Rear, the system can be in either magnetic inhibition regime 
Rm, or in the stage of supersonic propeller 
Rm (see, e.g., Bozzo et al. 2008). To de- 
cide between these latter possibilities it is useful to consider 
the relative extent of Reap with respect to i?/^- Indeed, when 

we akeady know that the matter proceeds within the magne- the pEase space parameters, and for a circular orbit (zero ec 

centricity). Several examples are marked (for different values 
of the magnetic field), represented by the solid lines. Each of 
these lines separate the behavior of the system from being a 
propeller (below it) to being magnetically inhibited (above it). 
We particularly note this division for B = 5 x lOi'^ G. Again, 



(1 -ecos(e)) ^) ^. (21) 



Fig. 10 shows the condition upon Voo for different values of 



tosphere, i.e., Rm < Ric, if Ric < Reap, we also have that 
Rm < Reap and the system behaves as a supersonic propeller. 
The inequality Rie < Reap happens when the period of the 
neutron star fulfills the following relation. 



Is/ I 108 cms-1 



lAMa 



1 - 0.69 



R* 



VlOi?c 



(l0i2 cm) 



the eccentricity value of the orbit is of importance. Fig. [11 
shows the ratio of the maximal Voo value below which tEe 
system is a propeller for a given eccentricity, with respect to 
the value it attains in a circular orbit. For instance, for e = 0.6 
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Fig. 10. — Phase space conditions for accretion regimes around a liighly 
magnetic neutron star. Above each of the curves, corresponding to different 
values of the magnetic field (an example is given for i? = 5 X 10^^ G) the 
system is a super-Keplerian magnetic inhibitor. Below this line, the system 
acts as a supersonic propeller. The eccentricity of the orbit is assumed as zero 
in this plot. The light-green shadow stands for the fact that the system halts 
matter within the magnetosphere in the whole of the phase space depicted. 




0.2 0.3 0.4 0.5 
Orbital eccentricity e 



0.7 



Fig. 11. — Ratio of the maximal Voo value below which the system is 
a propeller, for a given eccentricity, with respect to the needed value in an 
circular orbit. The higher the eccentricity the larger is the value of Voo that 
can accommodate a propeller behavior for a highly magnetic neutron stai'. 



the limiting velocity is enlarged by a factor 1 .65, meaning that 
each of the curves in Fig. [TOjshould be displaced in the y-axis 
by this factor. For the noted example with B = 5 x 10^^ G 
and the fiducial value of A/* — 10^^ g s^^, all systems with 
Voo — 2 X 10^ cm s^^ or less are supersonic propellers. 

4.3. Describing an equatorial disc outflow and caveats 

We will now consider the situation in which a neutron star 
feels the influence of an equatorial disc outflow, the matter 



in which is moving with velocity vector T4),eg- The uncer- 
tainties in describing the equatorial disc have a much more 
pronounced impact in the derivations. This section intends to 
give some of the details about this, and, if nothing else, make 
the caveats in the usual assumptions made in doing analytical 
treatments more explicit. 

4.3.1. A Bondi-Hoyle approach and its problems 

We start by coming back to Eq. (jTji for the Reap definition 
and first study which are the components that intervene in this 
scenario to give rise to the value of Vrei ■ Once Vrei is deter- 
mined, most of the (analytical) literature in the topic continue 
to assume that the accretion rate is given by 



Mace — TT 



{2GMnsf 



Peq = T^Rlap Vrei Peq, (22) 



where p^q is obtained from the continuity equation for matter 
radially outflowing, with velocity V^q^r, ^ disc of half-opening 
angle 9, 

M^" ^47rr^sm{e)Veq,rPeq, (23) 

and where M** is the mass loss rate in this outflow. These 
equations already imply that Bondi-Hoyle approximation is 
a valid one for the accretion process onto the star resulting 
from the equatorial disc outflow, i.e., that the material en- 
tering a cylinder of radius equal to the accretion radius will 
be captured. This may not be true by several reasons. First, 
because even when the Bondi-Hoyle could be assumed, the 
accretion radius can exceed the vertical size of the disc, and 
if so, the real accreted mass should be smaller (see, e.g., Za- 
manov 1995, Zamanov et al. 2001, and the discussion below). 
Also, short-lived Roche-Lobe overflows can occur in close bi- 
naries. This approximation also assumes that the Be disc is 
not affected in any way by the passage of the compact object, 
which is wrong based on simulations by Okazaki et al. (2002). 
In general. Be disks in binaries are tidally truncated (Okazaki 
& Negueruela 2001) and because of it, it is very hard to esti- 
mate the accretion rate onto the compact companion without 
running numerical simulations (see, e.g., Okazaki et al. 2002 
and Romero et al. 2007). 

It is nevertheless instructive to see how large an accretion 
rate one gets by using the Bondi-Hoyle approximation, and 
particularly, how much it is sensitive to changes in the as- 
sumptions made. Regarding the other components of Vrei 
apart those contained in V^.eq, from Kepler's law the orbital 
velocity of the neutron star around the star is given by. 



GM3 



(24) 



where is the mass of the stellar companion and the re- 
maining magnitudes have been already defined. In terms of 
the true anomaly 0, defined with the convention of having 
61 = at periastron, we can write 



a(l 



1 + 6008(6*)' 



(25) 



for the neutron star - Be star separation. As we will see, the 
magnitudes for the forthcoming assumptions for the equato- 
rial wind velocity make Vorb no longer negligible with respect 
to the other components of Vrei- We have already used that 
the equatorial disc matter will be assumed to move with a ra- 
dial, Veq^r velocity. In addition it may have also an azimuthal. 
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Veq.,p, component -assumed positive for rotation in the direc- 
tion of the orbital motion- such that V^ ,,^ — V^q.r + Kg </.■ 
Once these are defined, the relative velocity is given by t^ei — 
Vw,eq — Vorb ^nd its magnitude is given by 



* rel * orb ' * eq.,<p ' * eq,r 

-2VorbVeq,4, COs(0) - 2VorbVeq.r sin{(f>) , (26) 



where (f> is the flight-path angle between the neutron star ve- 
locity and the local horizon; measured from the latter to the 
neutron star velocity vector. The flight-path angle is positive 
when the neutron star is traveling from periastron to apastron, 
and negative otherwise. 

To proceed further, we need to define the functional form 
adopted for Veq^r and Veq.^- However, the magnitude of those 
is far from clear, as is the case, for instance, also for the ter- 
mination distance of the equatorial wind, or its granularity. 
In this exploratory Section, we can only base our study in 
different analytical formulations that have been already used 
in the literature to represent these components. The first ob- 
vious assumption is to consider that there is no azimuthal 
movement in the matter of the equatorial disc, Veq,<^. o = 
(e.g. Zamanov et al. 2001, Gregory and Neish 2002). Gre- 
gory and Neish (2002) considered several other parameteriza- 
tions for Veq,^, that go from a Keplerian velocity Vgq.^, Kcp = 
■\/ Gil/* / r; which would be the natural choice for a viscous 
decretion disc, to different power-laws assumed for the rota- 
tion of the envelope 14^,0. Men = V^,{r / R^,)^^^'^ (see also, 
Mennickent et al. 1994) where is the stellar rotation ve- 
locity -set to 360 km s^^, or Veq,0,pi = Vn^^{R^/r)/ sin(i) 
with i the inclination and Vb,^ ^ 1 km s^^ (see also Casares 
et al. 2005 and Bosch-Ramon et al. 2006). As can be seen 
in Fig. [12] these latter alternatives produce similar (within an 
order of magnitude) values for the magnitude of the azimuthal 
velocity component along the orbit, all of them also similar to 
the orbital speed. 

To define the radial velocity component one can also as- 
sume, to first approximation, that it is close to zero. Outflow 
velocities are so low that typically only upper limits (of the or- 
der of 1 km s^^) are obtained observationally (Okazaki 2010). 
Marlborough et al. (1997) used for instance Vo=0.3 km s^^ 
for their disc parameters. One could also adopt the power-law 



wind density p, 



eq 



with n = 3.2, used by Waters (1986) 
and Marti & Paredes (1995) in fitting their near-infrared data. 
For a constant outflow rate, the continuity equation leads to a 
radial outflow of the form Veq ,. — Vb(?'/i?s)"^ 2, where again 
r is the distance measured from the companion star, and Vo is 
quoted from a few to a few tens of km s^^. Vq was assumed 
as 5 km s^^ in Waters (1986) and in the range 2-20 km s"^ 
in Marti & Paredes (1995); what fixes the density of the equa- 
torial wind matter at the star surface on the order of 10^^^ g 
cm^'^, with a half-opening angle O 15°. Gregory & Neish 
(2002) obtain the radial velocity starting from the azimuthal 
one, forcing the accreted mass to be related to the radio emis- 
sion along the orbit. The radial velocities obtained in this 
manner are, for distances smaller than 10 i?s, quite similar to 
the functional form assumed above (see figure 5 in Gregory 
& Neish's paper). In addition, we take into account that the 
differences beyond this distance will not play a relevant role if 
the disc truncates, what is modeled with a phenomenological 
cut to the density of the equatorial wind beyond a distance of 
^ 10-R* (see, e.g., Gregory & Neish 2002, Bosch-Ramon et 
al. 2006). By assuming a cut at ~ 12i?* the second accretion 
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Fig. 12. — Components of the velocity vectors considered. In red (three-dot 
dashed) we plot alternatives for the assumption of the circular velocity of the 
wind matter In blue (dashed) we plot a couple of radial velocity components 
with different pre-factor. The orbital speed of the neutron star is depicted by 
the solid-black line. 



peak seen in some of the curves below is also preserved. 
In Fig. 



13 we plot Vrei as a function of the orbital phase 
-defined in terms of the eccentric anomaly as (/)p = (e — 
e sin(e))/ (27r), so that the phase is equal to zero at periastron- 
for the different alternative descriptions of the outflow previ- 
ously considered. In particular, we see that large differences 
appear along the orbital phases of the LS 1 +61° 303 system 
for the Keplerian's, Mennickent's, and power-law descrip- 
tions of the azimuthal velocity. For the latter case, this plot 
shows also how large is the influence of a changing initial ve- 
locity in Veq,r, ffom 3 to 5 km s^^ along the whole orbit. For 
completeness we also show the case of zero azimuthal veloc- 
ity (the case of zero radial velocity would produce a similar 
curve to the Keplerian case, depticted in red in the figure, par- 
ticularly at low phases). The accretion rate is very sensitive to 
small changes in Vrei- 

Fig. [13] also shows the accretion rate onto the neutron star, 
obtaine3~via Eq. (23 1, for a few of the cases. The solid lines 
stand for the results corresponding to the non-dissolving equa- 
torial disc; whereas the dashed lines correspond to those in 
which the disc density has been cut off beyond 10 Rs- The 
two blue lines correspond to the result with an initial radial 
velocity of 3 (top curve) and 5 km s^^ (bottom curve), re- 
spectively. The horizontal line stands for the stellar mass- 
loss rate assumed (1.3xlO~^ M0 yr~^), consistent with the 
equatorial disc density at the surface being (~ 1 x 10~^^ g 
cm^'^). The fact that some of the curves are on top of the 
horizontal line in Fig. [T3]must be understood as an unphysi- 
cal effect produced by tEe Bondi-Hoyle formula. This situa- 
tion was also found earlier (see, e.g., Zamanov 1995, Marti & 
Paredes 1995, Bosch-Ramon et al. 2006, Orellana & Romero 
2007, etc.) albeit the plotting of the ratio of the accretion rate 
with respect to the value it attains at periastron, might perhaps 
make this fact less clear at first sight. The failure of the Bondi- 
Hoyle formulae was treated differently by different authors, 
from ignoring it, given the approximate character of the stud- 
ies (see above-quoted references), to the addition of one or 
several ad-hoc cuts (Zamanov et al. 2001), to the use of radio 
data to re-scale it. In the last case, Gregory & Neish (2002) 
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assumed that the accretion rate value should be corrected by 
making S{r) — K{r)Macc, where S{r) is the radio flux den- 
sity (as measured along the orbit at e.g., 8.3 GHz) and K{r) is 
the factor of proportionality, depending on the objects separa- 
tion or orbital phase. The former expression can be re-written 
as k{r)p{r) = S{r)V^^i, where k{r) = K{r)'!i{2GMf . If 
Vrei is computed, S{r) is taken from data, and p{r) is as- 
sumed equal to Pf,q oc r^" as above, k{r) can be derived and 
K{r) immediately follows (Gregory & Neish 2002, see their 
figure 6). From this approach all values of the accretion rate 
are diminished significantly, and correspond to the range from 
0.001-0.01 Eddington's; which is equivalent to a reduction of 
>100 times the value of the accretion rate at periastron. This 
order of magnitude for the accretion rate would be consistent 
with the one numerically obtained by Romero et al. (2007). 
As noted by Orellana et al. (2007, see their figure 1), some 
of the accretion rate curves of Fig. [T3]display similar features 
-e.g., the position of the local maxima in phase- with the one 
numerically obtained, but there are approximately three or- 
ders of magnitude of difference in absolute value. We note 
that, consistently, Romanova et al. (2003) and Toropina et 
al. (2006) found that the fraction of the Bondi accretion rate 
which accretes to the surface of the star, obtained from MHD 
simulations of the accretion to a rotating star in the propeller 
regime, is much less than Bondi's. We come back to comment 
on these simulations below. 

4.3.2. Regimes and orbital eccentricity ejfects 
under the influence of an equatorial disc outflow 

The analysis in the previous section emphasizes that the 
problem of matter exchange in an equatorial wind setting is 
hard to quantify in an analytical treatment. It can likely be 
only assessed via numerical simulations, albeit these do not 
lack their own complexity and caveats (see, e.g., Zdziarski et 
al. 2010 and Okazaki et al. 2011 for an assessment). Never- 
theless, the discussion made in the last section is enough for 
the point we want to make here: at periastron distances and 
even considering a significant reduction of the accretion rate 
shown in Fig. 13 if the neutron star pass through the disc 



it is likely that accretion from the equatorial outflow dom- 
inates over that from the polar wind, implying M, 



acc,eg 



> 



M, 



acCjpolar- 



For comparison, we recall that t he p olar outflow 



accretion, computed by use of Eqs. (|6]) and ( 10 1, lead to val- 
ues that are maximal at periastron, where it attams ^ 2 x 10^"' 
g s^^ (this value is two orders of magnitude less than the ex- 
pected equatorial-outflow accreted-matter resulting from nu- 
merical simulations), to quickly fall to 5 x 10^^ g s^^. In this 
situation, it is even easier for a magnetar to allow matter en- 
tering into its magnetosphere, given that Rm is smaller under 
the equatorial wind influence. Indeed, 



i?„, = 8.2 X lO-" 



R \ 4/7 / ,V \ -2/7 
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-1/7 



12/7 



(27) 



and in order for matter in this outflow to enter within the mag- 
netosphere. 



1.4M, 







^1/7 



R„ 



106 cm 



12/7 



(28) 



The larger the accretion rate at periastron, the easier it is for 
pulsars with lower periods to have matter within the magneto- 
sphere, which is reflected, for fiducial values, comparing Eqs. 
([28) and ([B). 

Aditionally, to ask for the inequality _R„, > Rcor to be ful- 
filled implies, using Eqs. ( 18 i and (27 1, a condition over the 
period 
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which is easy to satisfy for all measured values of P of mag- 
netars. Finally, using the definition of Rcapi the condition 
Reap > Rm implies 
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(30) 



We lack accurate knowledge of Vrei to assess the former in- 
equality, other than observing that the fiducial value chosen 
is -to be conservative- one order of magnitude larger than 
the one shown in Fig. [13] A smaller value for Vrei (smaller 
than the one chosen asfiducial) would make the inequality 
even easier to fulfill. Then, this condition is not restrictive for 
magnetars, implying that, under the influence of an equato- 
rial outflow with an accretion rate of the order of 10^^^ g 
s^^, and with the caveat of not having a precise description of 
short-lived effects such as a Roche-Lobe overflow or a forma- 
tion of a transient accretion disc, the system would also act as 
propeller in its periastron. 

5. ENERGETICS 

5.1. Maximal energy of electrons at periastron shocks 

To analyze the possible maximal electron acceleration in 
the propeller we consider, as, e.g., Bednarek (2009), the rate 
of energy increase and of energy losses in a shock formed 
at the magnetic radius position; where turbulent motion of a 
strongly magnetized medium is considered prone to particle 
acceleration. The acceleration gain can be parameterized as 



dE 
~dt 



CcE/Rl = (ceB 



(31) 



where C, is an efficiency of acceleration (fiducial value to be 
taken as 10%), Rl is the Larmor's radius, and e is the electron 
charge. B is the magnetic field, which value will be obtained 
from Eq. (4]). ( isa free parameter here and the fiducial value 
for it has been chosen conservatively; generally <C 1, and 
is close to 1 only in extreme accelerators (e.g., Aharonian et 
al 2002, Khangulyan et al. 2007, 2008, considered values of 

C - 10^2 _ 10-4), 

Given that electrons of the highest energies loose energy 
mostly via the synchrotron process. 
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Fig. 13. — Left panel. Relative velocity along the LS I +61°303 orbit (phase zero is the system's periastron) for different alternatives for the description of 
the equatorial wind outflow. No density cut is assumed in this plot. The blue lines correspond to assume Vgq 0^ pi = 1.17 X 10^ ((iJ*/r) / sin(i)) cm s~^, for 
different assumptions of the Vo-factor in Veq.r- Righ t panel. Accretion rate onto the neutron star obtained by means of the Bondi-Hoyle formulation for different 
relative velocities. The color coding follows Fig.|13| see text for details. 



Here, ctt is the Thompson cross section, 7 is the electron 
Lorentz factor, and ps = B^/Svr is the energy density of the 
magnetic field. The maximum energy of the electrons are de- 
termined by the balance of the two former equations, which, 
by means of Eq. ([T4|, results in 
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where 



/(e) = (l-ecos(e))4/^ 
1-0.69^ ^* 



)(: 



IORqJ VIO12 cm 
(1 — ecos(e)) 1 



(34) 



In periastron, and with a — 6 x 10^^ cm, for fiducial values, 
the maximal electron energy prefactor (in Eq. (33 1) reaches 
6.8 X lO^'^ eV. If we instead use the equatorial result derived 
in Eq. ( 27 1, the maximal energy would read 



7(„^:?l-2.7xl0'^(^^ 
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(35) 



The maximal energy in this case is then of the order of 
1.3 X IQi^ eV, for fiducial values. Note that the higher the 
accretion rate (or the lower the magnetic field), the lower the 
maximal energy. Of course, the smaller the value of the effi- 
ciency, which could also change along the orbit, the smaller 
the maximal energy. To be explicit, for an assumed value of 
ry = 0.01 (e.g., Khangulyan et al. 2008), the maximal energy 
is reduced to 411 GeV. The turbulent region of the propeller 
shock would not make electrons to achieve TeV energies. 

6. DISCUSSION 

In the previous Sections we have discussed the conse- 
quences of having a magnetar in the LS 1 +61° 303 system 
(or at least a high B, long P pulsar able to burst in the way 
described). This possibility was based: 

• on the observational evidence for one very short flare 
observed by Swift, for which the only remaining coun- 
terpart candidate X-ray source within an improved Icr 
uncertainty in position is LS I +61° 303, 

• on the complete similarity of the previous event with all 
properties known for other magnetar flares, 

• and on the fact that the longer flares observed 
by other satellites (and undoubtfully coming from 
LS 1 +61°303, having timescales of about 1 h) are un- 
correlated with harder X-ray emission, so discarding 
that the shorter flare is a result of a strong spectral evo- 
lution of a typical non-magnetar burst (proven by an 
uncorrected 5w//f-BAT or /JXr/s-HEXTE detection in 
any of the longer flares). 

If we are to accept that a magnetar may be part of the 
LS I +61° 303 system (which would then be the first such 
in being discovered) we concluded from the analysis in the 
former Section that it would likely be subject to a flip-flop be- 
havior. It would shift from behaving as a neutron star being 
rotationally powered near apastron to being a propeller near 
periastron, along each of the system's orbit. 

Based on the preliminary analysis of the ^wZ/f-BAT burst as 
reported in Astronomer Telegrams (De Pasquale et al. 2008; 
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Barthelmy et al. 2008; Dubus & Giebles 2008), earlier con- 
siderations ofLSl+61°303being formed by a magnetar have 
been made (Bednarek 2009, Dubus 2010). In them, the sys- 
tem has not been proposed to flip-flop on any behavior, and 
the influence of the orbit was not considered. Dubus (2010) 
proposed that the system would always behave as a normal 
pulsar, mimicking in all respects the properties needed by an 
inter-wind shock to be sustained (e.g., as in Dubus 2006). He 
also proposed that the short burst suggested a magnetic field 
strength of 10^'' G, which we do not find justified. Instead, 
Bednarek (2009) proposed that the LS 1 +61° 303 system is al- 
ways a propeller, even in apastron. We find above that for typ- 
ical magnetar parameters, neither of these possibilities seems 
likely to be the case. 

The idea of a flip-flop behavior for a neutron star in an ec- 
centric orbit was earlier considered (see, e.g., lllarionov & 
Sunyaev 1975, Gnusareva & Lipunov 1985, Lipunov 1987, 
Lipunov et al. 1994, Campana et al. 1995). And actually, a 
flip-flop behavior was even proposed for LS 1 +61° 303 itself 
(Zamanov 1995, Zamanov et al. 2001) although the idea was 
quickly abandoned in the literature perhaps because for ms- 
pulsars, the system would indeed always be rotationally pow- 
ered, as shown above. In this work, we have thus rekindled 
some of these ideas but given them an extra edge: prompted 
by the observational analysis, the LS 1 +61° 303 pulsar has 
magnetar parameters and thus the flip-flop behavior is rather 
an expected outcome of the orbital evolution. But may the 
flip-flop behavior be in qualitative agreement with the remain- 
ing multi-wavelength phenomenology of the system? 

6. 1 . A flip-flopping magnetar model for LS I +61° 303 in the 
context of multi-wavelength observations 

6.2. TeV range 

The TeV emission of the LS 1 +61° 303 system has been 
discovered by MAGIC (Albert et al. 2006). A claim of pe- 
riodic recurrence of the TeV emission was made thereafter, 
with the system showing regular outbursts at TeV energies in 
a broad-range of phases around 0.65 (radio phases, as in Gre- 
gory 2002, where periastron is at phases 0.23-0.3, see, e.g., 
Casares et al. 2005 or Grundstrom et al. 2007), with no sig- 
nificant signal elsewhere (Albert et al. 2009). The VERITAS 
array did independent observations and soon confirmed the 
MAGIC detection (Acciari et al 2008, 2009). In addition, the 
MAGIC collaboration claimed a correlated X-ray / TeV emis- 
sion at the position of the outburst (Anderhub et al. 2009) 
which would emphasize the likely common origin of the ra- 
diation in both bands. However, this was based on 60% cov- 
erage of just one orbit, and there was no clue on how general 
this result really was. 

Neutron star-based models presented to analyze the TeV 
emission from LS 1 +61°303, would produce a recursive or- 
bital behavior, unless of course significant parameters (e.g., 
those defining the level of absorption) change in an orbit- 
to-orbit basis for yet unknown reasons (see, e.g., Dubus 
2006b, 2010; Sierpowska-Bartosik & Torres 2007, 2008, 
2009; Bednarek 2009). The strong dependence on the sys- 
tem's geometry and the orbital position shown by the absorp- 
tion and emission process (pair production and anisotropic in- 
verse Compton scattering) would produce an orbit-to-orbit re- 
current TeV burst. However, the peak of the TeV emission 
found by early measurements of MAGIC and VERITAS, at 
radio-phase 0.6 — 0.7, is not strictly near the superior con- 
junction of the orbit, where the conditions would be more fa- 



vorable to gamma-ray production. Additional considerations 
where made to correct for this displacement in these mod- 
els. Dubus et al. (2010) consider the influence of a relatively 
mild Doppler-boosted emission; assuming that there is an out- 
flow velocity of 0.15-0.33 c, consistent with the expected 
flow speed at the pulsar wind termination shock. Sierpowska- 
Bartosik & Torres (2009) assumed that the lightcurve is gen- 
erated by a variable power in the injection of relativistic par- 
ticles. 

Nevertheless, new observations by VERITAS (Acciari et al. 
2010) discard the regularity of the TeV burst and the X-ray / 
TeV correlation. The observing seasons from 2008 to 2010 
covered a total of 8 orbits out of 27 available. In none of them 
the system was detected where it was earlier expected, near 
apastron. Actually, in 7 out of 8 orbits the system was not 
detected in any phase; which can not be accommodated, as 
discussed by Acciari et al. (2010), arguing for shorter obser- 
vation times. The upper limit imposed were less than 5% of 
the Crab Nebula flux in the same energy range, in contrast to 
previous observations by both MAGIC and VERITAS which 
detected the source during these phases at 16% of the Crab 
Nebula flux. This implies that the system is in a low TeV state 
where the flux is at least a factor of 3 lower than the one de- 
tected in the 2006 observing season. Additionally, VERITAS 
reported the system was significantly detected just once, dur- 
ing observations taken in late 2010, at a phase much closer to 
periastron passage. 

In the scenario of this paper, the usual appearance of the 
TeV emission near apastron would be explained under the 
common lore: the inter-wind shock produced by the interac- 
tion of a rotationally powered pulsar and the stellar wind of 
the star accelerates particles, which in turn emit gamma-rays. 
The common disappearance of the TeV emission in perias- 
tron, instead, is accommodated by the fact that on the one 
hand, the energy cutoff for particle acceleration in the shock 
formed in the propeller regime is sub-TeV (and will be less for 
larger accretion rate, smaller magnetic field, and the expect- 
edly smaller acceleration efficiency out of the fiducial values 
assumed); and, on the other, the cross section for pair produc- 
tion for any TeV particle is maximal. Typically, then, the sys- 
tem would appear and disappear where it has been observed 
along the 2006 to 2009 seasons. 

Contrary to other models, the flip-flopping LS 1 +61° 303 
system could also provide a natural interpretation for the flux 
reduction or the total disappearance of the TeV radiation. If 
the mass-loss rate (and thus the accretion rate onto the neu- 
tron star) increases, the system more quickly abandons its 
rotationally powered behavior, and if it increases enough, it 
would not be in such a behavior at any portion of the orbit. 
The inter-wind shock would not form and abundant TeV par- 
ticles would not be produced. We note -see Fig. [9]- that even 
a small factor increase in the mass-loss rate at apastron can 
make the system to have a disrupted magnetosphere there as 
well. (If the mass-loss rate of the equatorial outflow increase, 
its properties may change significantly, e.g., the truncation ra- 
dius could be enlarged and even dominate at apastron, also 
implying that significantly less TeV emission would be pro- 
duced there.) 

The orbit-to-orbit TeV changes can be understood in equal 
terms, by directly linking it to the variations of the mass-loss 
rate in a highly variable behavior of Be stars. This has the ob- 
vious caveat of making a particular orbit's TeV lightcurve im- 
possible to foresee, albeit links observations able to constrain 
the mass-loss rate with those at the highest energies, and al- 
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lows for a prediction of anti-correlation to be made: the higher 
the mass-loss rate the lower the TeV emission. Allowing one- 
self to speculate, detailed observations of this anti-correlation, 
if confirmed, can impose constraints over the magnetar pe- 
riod. 

6.3. The TeV emission and its apparent correlation with the 
radio-obtained super-orbital phase 

We now compare the years-long evolution of the system 
as tracked in radio and TeV energies and look at it from the 
perspective of the magnetar-composed LS 1 +61° 303 model. 

For the radio data we use the measurements compiled by 
Gregory (1999), most of it at a frequency of 8.3 GHz, us- 
ing the Green Bank Interferometer The analysis of the radio 
emission led to the discovery of a super-orbital period at 1667 
days (Gregory 2002). Fig. [14] shows these data folded with 
the super-orbital period. Super-imposed, we also show all 
the super-orbital phase-ranges where there were TeV obser- 
vations (either by MAGIC or by VERITAS) which covered the 
orbital phases where the LS 1 +61° 303 system was initially 
detected near apastron, i.e., at a broad radio phase around 
^ 0.65. Then, all colored boxes in Fig. [T4| represent the 
TeV observations that covered the broadly-oeiined apastron 
region of the LS 1 +61° 303 orbit. Typically, several orbits 
are included in each of the colored boxes. In order to plot 



these TeV observations in Fig. 14 the MJD of each observa- 
tion time that covered the apastron region was obtained from 
all TeV reports on LS 1 +61°303 (Albert et al. 2006, 2009; 
Anderhub et al. 2009; Acciari et al. 2008, 2010) and con- 
verted into super-orbital phase using the period of 1667 days. 
We note that, although sparse, the TeV coverage of the source 
goes from 2006 to 2010. 

From Fig. 14 we see that the system has appeared as a sig- 
nificant TeV-source near apastron only (light yellow boxes) at 
super-orbital phases near the minimum of the radio flux. In- 
stead, at the maximum of the radio flux, i.e., at super-orbital 
phases near ^ 0, observations confirmed only TeV upper lim- 
its (green colored boxes). The latter situation corresponds to 
the recent disappearance of the source as reported by VERI- 
TAS, signifying an important flux reduction. Needless to say, 
the radio data and the TeV measurements are not contempo- 
raneous, and we will especulate that the average properties of 
the super-orbital period producing the radio long-term modu- 
lation have been maintained since 2002 till now to extract con- 
clusions. This may not necessarily be true (see, e.g., Trushkin 
& Nizhelskij 2010), although the drift in the maxima of the 
radio peak seems very small (of the order of 3 days) to affect 
the reasoning. 

In itself, it is possible that Fig. 
the nature of the source. If we follow the interpretation by 
e.g., Gregory et al. (1989), Zamanov et al (1999), and Gre- 
gory (2002) and accept that the super-orbital radio flux mod- 
ulation is the result of a cyclical variation in the mass-loss 
of the Be star (the synchrotron radio power emitted is di- 
rectly proportional to the relativistic particle density which 
likely scales with the mass-loss rate), the disappearance of 
flux-reduction of the TeV radiation at the peak of the radio 
emission can be accommodated as follows. In the magnetar 
model for LS I +61°303, the increase in mass-loss rate leads 
to an increase in the accretion rate onto the compact object, 
which pressure makes the magnetic radius to become similar 
and eventually less than the light cylinder, disrupting the mag- 
netosphere even at apastron. If so, the system stops behaving 
as a normal pulsar, stops driving a wind, and stops forming 
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an inter- wind shock. Instead, LS I +61° 303 acts as a pro- 
peller. As a result, it is not expected to produce a significant 
number of multi-TeV particles, and thus the TeV radiation is 
suppressed. As in periastron the system was already behaving 
as a propeller, there is no significant TeV in any part of the 
orbit. 



We note that the results of Fig. 14 are model-independent 



since they are based only on observational results (when there 
was a detection near apastron, when there was none) and the 
times in which these observations were taken. They are how- 
ever difficult to accommodate in other models of the source 
(which in addition are not able to explain the short burst). On 
the one hand, in a microquasar scenario, the increase of the 
mass-loss rate would enhance -under the disc-jet coupling as- 
sumption (e.g., Falcke & Biermann 1995)- the power in the 
jet. Thus, in a microquasar model, when the mass-loss and 
accretion rate increase there should be more TeV radiation, 
not less. On the other hand, estimations of the variations in 
the mass loss-rate from the Be star are given as the ratio be- 
tween maximal and minimal values obtained either from radio 
emission (a factor of 4 was determined by Gregory & Neish 
2002) or from Ha measurements, which span from a factor 
of 5.6 (Gregory et al. 1989) to 1.5 (Zamanov et al. 1999); 
in any case, a factor of a few. In a normal pulsar wind - 
stellar wind scenario, when the pulsar is not a magnetar, an 
increase in the mass loss rate by a factor of a few does not 
change anything in the behavior of the source (what can be 
verified using for instance Eq. ([T5|)): the system would al- 
ways have a non-disrupted magnetosphere. In the inter-wind 
shock scenario with a normal pulsar being part of the system, 
there is no apparent reason for the disappearance or the sig- 
nificant reduction of the TeV emission, nor for the apparent 
correlation of this reduction with the maximum of the mass- 
loss rate when it increases by a factor of a few. Instead, in 
the case of a magnetar in LS I +61°303, fiducial values of the 
involved magnitudes put the neutron star in a position of the 
phase-space in which small changes in the mass-loss rate can 
produce regime changes, accommodating Fig. 14 



6.4. GeV range 

GeV emission from LS I +61° 303 has been detected using 
the Large Area Telescope onboard Fermi (Abdo et al. 2009), 
at the orbital periodicity. With the data of the first 8 months of 
the mission, only a single spectrum averaged along the orbit 
could be obtained. A power-law with an exponential cutoff 
at a few GeV was found as the best fit. The folded Fermi 
lightcurve peaked at periastron passage. The anti-correlation 
of the phase of the maximum between TeV and GeV can be 
understood as a result of inverse Compton scattering and pair 
absorption. With 2.5-years of Ferm/-LAT data more details 
have been obtained (see, e.g., Torres et al. 2010, Hadasch 
et al. 2010 for preliminary reports). The spectral characteri- 
zation of LS I +61° 303 allows now for a fit of a power-law 
with an exponential cutoff spectrum along each analyzed por- 
tion of the systems orbit (which is typically cut in halves). 
The differences between these spectra is however not large, 
of the order of 20% in flux and within errors in the determina- 
tion of the cutoffs. The LAT now detects emission up to tens 
of GeV, where prior datasets led to upper limits only. Due 
to the contemporaneous measurements by VERITAS, com- 
mented above, we know now that there is no evidence that 
the process responsible for the detected Fermi-LAI emission 
continues beyond the GeV cutoff. One of the most interesting 
results of this further Fermi-LAT campaign is the detection of 
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Fig. 14.— Radio data from LS I +61°303 compiled by Gregory (1999), 
most of it at a frequency of 8.3 GHz, and folded with the super-orbital period 
of 1 667 days. The colored boxes represent MAGIC and VERITAS observation 
times that covered the broadly-defined apastron region of the LS 1 +61° 303 
orbit. Light-yellow boxes stand for periods in which the LSI+61°303 sys- 
tem was detected whereas the light-green boxes stand for those in which it 
was not, implying at least a factor of 3 reduction in flux. Each of the boxes 
represent several orbits of LS 1 +6 1 ° 303 . Two super-orbital cycles are shown 
for clarity. 



an overall flux increase around March 2009, of the order of 
30%. This flux change was accompanied by a significant de- 
crease in the GeV lightcurve modulation, which since then is 
much flatten. This flux change is then approximately coinci- 
dent with the period of TeV flux reduction. 

Within the magnetar-composed flip-flopping model of 
LS 1 +61°303 the usual behavior of GeV emission, being 
modulated along the orbit and anti-correlated with the TeV 
emission, finds its place as commented before. The increase 
in GeV luminosity like the one found after March 2009 can 
be understood due to the increase of the accretion rate, which 
would in turn increase the energy reservoir of which a fraction 
ends up in the GeV domain. If this is the case, the magne- 
tar in LS 1 +61° 303 would be behaving as a propeller along 
more portions of the orbit, and eventually along all of it. One 
would qualitatively expect that the GeV emission would not 
be significantly modulated anymore, since it is the same pro- 
cess generating the radiation all along the orbit (no more inter- 
wind shock), the GeV modulated fraction would diminish and 
flatter, with the TeV emission being simultaneously reduced. 

6.5. X-ray range 

Zhang et al. (2010), Torres et al. (2010) and Li et al. (201 1) 
present the status of our knowledge of the X-ray behavior of 
the source, via long-term, years-long monitoring using the IN- 
TEGRAL and RXTE satellites. The latter is the richest dataset, 
as has been discussed and enlarged via the addition of several 
more orbits in Section[T]above. Flares in the ks-timescale have 
been detected since long by different experiments (see refer- 
ences in the quoted papers) and they can be accommodated as 
local processes, for instance, through dumpiness in the stel- 
lar wind. Their interpretation would likely not be affected 
within the magnetar-composed LS 1 +61° 303 model. An 
outcome of the RXTE campaign is the discovery of variabil- 



ity in the orbital profiles of the X-ray emission. Lightcurve 
profile variability is seen from orbit-to-orbit all the way up 
to multi-year timescales, with the phase of the profile maxi- 
mum also changing. Such a high degree of variability can cer- 
tainly be more easily accommodated with a magnetar hosted 
by LS 1 +61°303 than in other models. In the former, we 
have different processes at different places contributing to the 
X-ray emission, already in an orbit-to-orbit basis, the predom- 
inance of which and the phases at which they become dom- 
inant can be altered by both, random (in orbital timescales) 
and cyclical (in super-orbital timescales) variations of the ac- 
cretion rate. The possible appearance of the super-orbital vari- 
ability in X-rays remains a matter of study for when additional 
data accumulates. 

6.6. Radio range 

The existence of periodic (^ 26.5 days) non-thermal radio 
outbursts have allowed to determine a strong modulation in 
the amplitude occurring on a timescale of ^ 4 years, both of 
which have been commented above (see, e g., Gregory 2002, 
and references therein). Cyclic variability in the Be star enve- 
lope have been proposed as a possible origin of the long term 
modulation. This has been emphasized by an apparent cor- 
relation of the super-orbitally varying radio flux with the Ha 
line parameters; the latter also present super-orbital modula- 
tion with the same period, albeit it is shifted in phased about 
400 days (e.g., Zamanov & Marti 2000). In terms of the flip- 
flopping magnetar model for LS 1 +61° 303, the outburst on- 
set could represent the time of the transition of the neutron 
star from propeller to a normal pulsar behavior (see, e.g, Za- 
manov 1995). The discussion in Zamanov & Marti (2000) put 
the general features of this model in context with radio obser- 
vations (note that the super-orbital period has been coiTected 
from 1584 to 1667 days since that paper was published, and 
the analysis should be changed accordingly). 

Very Long Baseline Array (VLBA) imaging obtained by 
Dhawan et al. (2006) over a full orbit of LS 1 +61° 303 has 
shown the radio emission to come from angular scales smaller 
than about 7 micro-arcsec (which is a projected size of 14 
AU at an assumed distance of 2 kpc). This radio emission 
appeared cometary-like, and interpreted to be pointing away 
from the high mass star and thus being the smoking gun of a 
pulsar wind. These results were found to be consistent with 
observations by the MAGIC collaboration, simultaneously us- 
ing MERLIN in the UK, the European VLBl Network (EVN), 
and the VLBA in the USA (Albert et al. 2008). The compari- 
son between Dhawan et al. (2006) and Albert et al. (2008) im- 
ages at the same orbital phase (but obtained 10 orbital cycles 
apart) show a high degree of similarity on both its morphology 
and flux, which suggests a certain stability of the physical pro- 
cesses involved in the radio emission. The tail is not always 
seemingly pointing in the right direction: as pointed for in- 
stance by Romero et al. 2007 or Zdiarszki et al. 2010, a high 
spin-down power pulsar would generate a wind which should 
overcome the stellar wind power, generating a flat inter-wind 
shock, or even wrapping it around the Be star, not the pulsar. 
Morphology-wise, the flip-flopping model may alleviate this 
problem, since even for a high E pulsar, when a propeller is 
running, there is no pulsar wind and the high-energy emis- 
sion is not coming from the inter-wind shock. The analysis 
of the conical shape of the radio emission would not apply. 
MHD simulations of accretion onto a magnetized neutron star 
in the propeller regime have been recently been presented by 
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Romanova et al. (2003) and Toropina et al. (2006). Their 
simulations showed that accreting matter is expelled from the 
equatorial region of the magnetosphere moving away from the 
star in a supersonic, disk-shaped outflow. Whether this out- 
flow may be related with the radio morphology is yet an open 
question. 

6.7. Absence of pulsations at all frequencies 

Finally, we remark that it is natural to expect a failure in de- 
tecting pulsations from a magnetar-composed LS 1 +61° 303 
system. Radio pulsations are not common in magnetars. 
A few examples of radio pulsed emission were detected 
mostly connected with the outbursts of a few transient ob- 
jects (Camilo et al. 2006, 2007; Levin et al. 2010), but it is 
far from being a common or stable property of these highly 
magnetized neutron stars. 

On the other hand, apastron would be the only region of 
the orbit when, for values of accretion rate so permitting, the 
pulsar in LS I +61°303 would have an unscathed magneto- 
sphere. However, the stellar winds might have prevented ra- 
dio pulsations to be detected because of the strong free-free 
absorption all over the orbit and/or the large and highly vari- 
able dispersion measure (DM) induced by the wind. This 
would be very similar to the case of PSR B1259-63, which 
at an inter-stellar distance that has about the same dimension 
of the major axis of the orbit of LS I +61° 303 (given the 
larger 3.4-yr orbit of that system), does not show radio pulsa- 
tions at periastron either. In X-rays, the upper limits derived 
by Rea et al. (2010) using Chandra observations at phases 
close to apastron, mimic the situation of PSRB1259-63 (see 
the discussion in Rea et al. 2011), where again no pulsation 
is found. Furthermore, X-ray pulsations from magnetars can 
have pulsed fractions as low as a 4% (e.g., Mereghetti et al. 
2007), well below the ^ 10% upper limit derived by Rea et 
al. (2010). The absence of X-ray pulsations can be also un- 
derstood in terms of having an X-ray emission which may 
be dominated by the shock rather than being due to emission 
from the pulsar magnetosphere. In periastron, and in general, 
in all phases where matter enters within the light cylinder, the 
magnetosphere is disrupted and pulsating radio emission is 
halted. X-ray pulsating emission, on the other hand, is also 
halted in this case because the accreting matter does not reach 
the surface of the neutron star. 

7. PREDICTIONS & OUTLOOK 

We can summarize a few testable predictions that are inher- 
ent to the scenario we have presented. Some of them are pre- 
cise, and some represent trends on which further study could 
shed light on the validity of the assumptions made: 

• If one can track the mass-loss rate evolution of the com- 
panion star in LS 1 +61° 303, the TeV emission will 
be anti-correlated with it. This would be valid at all 
timescales, e.g., both in an orbit-to-orbit basis, as well 
as in longer super-orbital timescales. 

• Assuming the persistence in time of the super-orbital 
radio modulation, if the cyclical increase of the mass- 
loss rate of the Be star already produce a flip-flopping 



LS 1 +61°303 system as suggested by Fig. 14 it would 
be natural to expect low TeV fluxes near apastron for 
super-orbital phases of ^ 1 ± 0.2. 

• Thus, within this model and under the assumption of 
the maintenance of the cyclical behavior of the accre- 



tion rate as inferred by the (somewhat old) radio data, 
it is predicted that the system should be already visible 
by TeV instruments. It should be back in the high state 
since approximately 2010 May - June, when the system 
attained super-orbital phase of 0.2 (there is no coverage 
of the apastron orbital phases at this epoch reported in 
Acciari et al. 2010) and will disappear or have the TeV 
emission severely lowered again at super-orbital phase 
0.7 - 0.8, or about 2012 October 

• If there is ever a giant flare observed from the magnetar 
hosted in LS I +61°303, it will aflow the detection of a 
magnetar-like pulsating tail, which will give us a handle 
on the rotational period of the neutron star 

• It would be natural to expect that when the TeV emis- 
sion is back to normal (appearing near apastron) the 
GeV to TeV connection there should be different from 
that in periastron (this could not be tested yet, since the 
TeV emission was reduced only 8 months after Fermi 
launch, and an orbital-separated spectrum could not 
be derived with such dataset). In apastron, the TeV- 
emitting electrons are present, which may lead to the 
appearance of a second component above 10 GeV in 
Fermi data. This should be less visible in periastron, 
where the maximal energy of electrons is cut at lower 
energies by the synchrotron losses. 

• It would be natural to expect that when the TeV emis- 
sion is significantly reduced or completely disappears 
along the orbit, including near apastron, the radio mor- 
phology should change, given that there is no pulsar 
wind in place and it is instead the propeller what drives 
it. 

We emphasize that a real description of these processes is 
surely far from the simplistic approach presented here. The 
transitions between the regimes are in fact not modeled, and 
will likely depend on how matter is accreted, on whether a 
temporary accretion disc is formed, and on many other details. 
These will probably be only resolved through numerical sim- 
ulations. This paper should rather be taken as an exploratory 
work for the plausibility of the idea when confronted with the 
rich phenomenology of the system. Under these caveats, we 
have found that the flip-flop magnetar LS I +61° 303 model 
not only seems to be suggested by a close look at the short 
burst observed by Swift-BAI and in general, at the X-ray 
variability, but appears able to qualitative agree with all other 
multi-frequency observations of the system and may explain 
the long term behavior of the TeV variability. This qualitative 
ability of the model to accommodate a wide variety of con- 
straints, even where others seem to fail, encourages further 
study. In particular, studies of accretion onto strongly magne- 
tized systems may prove useful. 
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